Vortex Glass and Vortex Liquid in Oscillatory Media 
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We study the disordered, multi-spiral solutions of two-dimensional homogeneous oscillatory media 
for parameter values at which the single spiral/vortex solution is fully stable. In the framework of the 
complex Ginzburg-Landau (CGLE) equation, we show that these states, heretofore believed to be 
static, actually evolve on ultra-slow timescales. This is achieved via a reduction of the CGLE to the 
evolution of the sole vortex position and phase coordinates. This true defect-mediated turbulence 
occurs in two distinct phases, a vortex liquid characterized by normal diffusion of individual spirals, 
and a slowly relaxing, intermittent, "vortex glass". 



PACS: 45.70.-n, 45.70.Ht, 45.70.Qj 

Spiral waves are ubiquitous in oscillatory and excitable 
two-dimensional active media Q . Their cores are robust 
wave sources which determine the oscillating frequency 
of the entire system and may dominate the surrounding 
dynamics. Spirals (often called vortices) may sponta- 
neously appear and annihilate in a typical manifestation 
of spatiotemporal chaos. When the single spiral solution 
is stable, one easily observes locked, quasi-frozen, multi- 
spiral disordered structures whose glassy character has 
been suggested ^|-^| . Despite being found in most models 
of excitable or oscillatory media as well as in experiments 
(Belousov-Zhabotinskii reaction H, surface growth ||, 
and others) surprisingly very little is known about the 
properties of these disordered states and the transitions 
leading to their formation. 

In this Letter we investigate, in the framework of the 
complex Ginzburg-Landau equation (CGLE), the general 
issue of these multi-spiral solutions and show that they 
actually evolve on ultra-slow timescales. This is achieved 
thanks both to long numerical simulations of the CGLE 
and to a quantitatively correct reduction of the dynamics 
to the evolution of the sole vortex position and phase co- 
ordinates. This true defect-mediated turbulence occurs 
in two distinct phases, a vortex liquid characterized by 
normal diffusion of individual spirals, and a slowly relax- 
ing, intermittent, "vortex glass" . 

The CGLE describes most properties of generic oscil- 
latory media, at least at a qualitative level, even if one 
is not in the vicinity of a supercritical long-wavelength 
Hopf bifurcation, where it can be systematically derived 
(for reviews, see fi]||). Under appropriate scaling of the 
physical variables, it takes the universal form 



d t A = A - (1 + ic)\A\ 2 A + (1 + ib)AA, 



(1) 



a wide variety of striking dynamical phenomena in one, 
two, and three space dimensions, many of which were also 
observed in various experimental contexts, sometimes up 
to a quantitative agreement with CGLE predictions. 

A distinctive feature of the two-dimensional CGLE is 
the existence of nontrivial sources of spiral waves (vor- 
tices) which determine the oscillating frequency of the 
entire system H. The single-spiral solution for a vortex 
centered at r reads as: pjjLl|: 



(2) 



A s (r, t) = F(r) exp[i(md + tp(r) + cut + <p)}, 



where A is a complex amplitude, b and c are real pa- 
rameters characterizing relative dispersion and nonlinear 
frequency shift, and A is the Laplace operator. Inten- 
sive studies conducted over the last ten years revealed 



where r — |r — ro|, is the polar angle measured from 
the vortex core, <fi is an arbitrary phase, and m = ±1 
is the topological charge. Far away from the core, the 
solution approaches a plane wave with ip(r) ~ kr, where 
the asymptotic wavenumber k is related to the rotation 
frequency as uj = —c — (b — c)k 2 . The dependence of 
k on b and c is known analytically for \b — c| -C 1 and 
\b - c\ > 1, e.g. k ~ -e _1 exp(— 7r/|2c|) for b = and 
\c\ -C 1 Jl(|. In the parameter region where the single- 
spiral solution is stable and c ^ 6, the interaction between 
two well-separated spirals falls off exponentially [|ll]-|l3j. 
The screening was attributed to the shock lines where 
the waves emitted by the cores collide. For relatively 
small b, c the interaction is monotonic, and the spirals ex- 
hibit asymptotic repulsion irrespectively of their charge 
("monotonic range") jll 12|. In contrast, for larger b, c 
satisfying the condition (c - b)/{\ + be) > c* ps 0.845 
the velocity vs distance dependence is modulated ("os- 
cillatory range" ) , and the spirals become keen to form a 
variety of long-living (but unstable) bound states |n[[L4| . 

While the above interactions cannot account for the 
strongly chaotic regimes where many defects are spon- 
taneously generated and undergo violent motion, they 
are expected to play a leading role in the occurrence 
of the quasi-static structures of large and small spirals 
surrounded by a complex network of shocks commonly 
observed in the large region of the (b, c) parameter space 
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where the single spiral solution is stable |2j-|4j . Rather lit- 
tle is known about these cellular structures except that 
they seem to be frozen on usual observation timescales. 
Our results lead to a dramatically different picture. 

We conducted detailed numerical studies of the two- 
dimensional CGLE restricted to the case 6 = and 
c > (for the regimes of interest scaling relations ap- 
ply in parameter space ||). The integration domain 
with periodic boundary conditions was typically of area 
S = 512 x 512 and the integration time about 10 . Vari- 
ous integration schemes were used and all results checked 
against variation of the numerical resolution. We mon- 
itored the positions Tj(t) of all N spirals. We mea- 
sured their instant "activity" T = 5 _1 J s dr \dt\A\\ and 
— when applicable — the spiral diffusion coefficient D = 
-^j 5Zj=i (l r i(*) — r j(0)| 2 )- Since, for well-separated spi- 
rals, \A\ is constant everywhere except near the cores, T 
is related to the velocity of spirals. 
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FIG. 1. Vortex liquid in the monotonic range, for the full 
CGLE (a,b,d,e) and for Eqs.(|^) (c,f). (a,b) snapshots of 
Ke(A) and \A\ for a typical solution of the CGLE at c = 0.6 
with N = 82 spirals (domain of size 512 x 512 with peri- 
odic boundary conditions), (c) Snapshot of core positions 
and shock lines for a typical solution of Eqs.(|^) for the same 
parameters as in (a,b) (d) Mean square displacement (r 2 ) 
of vortices vs time obtained from solutions of the CGLE for 
c = 0.7, L = 256, and N = 28,22,18 vortices from top to 
bottom, (e) Distribution function of the instantaneous spiral 
core velocity (arbitrary units) with fit by a Maxwellian from 
a simulation of the full CGLE at c = 0.7 with 118 vortices in 
a 512 x 512 domain, (f) Variation of D, the diffusion constant 
of vortices, with l/^/p for c = 0.3 (from solutions of Eqs. ([])). 



In the monotonic range c < c* (Fig. [|), after a short 
transient during which the initial number of defects may 
decrease, the number of spirals remains constant and T 
fluctuates around a well-defined mean value. The typi- 
cal vortex velocity is very small (of the order of 10~ 4 for 
c = 0.7), and the velocity distribution is approximately 



Maxwellian (Fig. Qd). In large enough systems, the spi- 
ral cores perform normal diffusive motion (Fig. |l|d). This 
vortex liquid can be characterized by a viscosity v ~ D^ 1 . 
Since all scales, already very large near c*, diverge expo- 
nentially when c — * 0, an extensive investigation of the 
dependence of D on c and on the vortex density p = N/ S 
is beyond reach of present-day computers. Nevertheless, 
our data indicate that D increases with p. 




FIG. 2. Results in the oscillatory range (c = 1.2), obtained 
on the full CGLE (a,b,d,e) and from the reduced equations 
(c,f). (a,b) snapshots of Ke(A) and \A\ for a typical solution 
of the CGLE with N = 124 spirals (domain of size 256 x 256 
with periodic boundary conditions), (c) Solution of Eqs.(^) 
for the same parameters as in (a,b). (d,e): time series of the 
activity T and of the number of vortices N for the evolution of 
the CGLE in a domain of size L — 256 starting from random 
initial conditions, (f): same as (d), but from simulations of 
Eqs.(|) with TV = 100 vortices (in this case T = N' 1 YV |vj|) 



In contrast, in the oscillatory range c > c*, after 
a very long, possibly infinite, transient (Fig. ^Js), the 
population of spirals spontaneously segregates into two 
distinct phases: large and almost immobile spirals and 
droplets (clusters) of small vortices confined between 
them (Fig. ||a,b, see also [^|-Q). One can thus define 
a "liquid fraction" of small spirals, whose sizes are typi- 
cal of the vortex liquid observed in the monotonic range. 
When the liquid fraction is small, the resulting state ex- 
hibits slow intermittent dynamics (bursts of activity sep- 
arated by long quiescent intervals, Fig. ||d) reminiscent of 
glassy dynamics, whereas for higher fractions the activ- 
ity T fluctuates around some mean value and the "liquid" 
vortices can be shown to perform normal diffusion in the 
space surrounding the big spirals. This picture is remi- 
niscent of phase coexistence in first-order transitions of 
isolated equilibrium systems Jlq| . In this respect, our re- 
sults indicate that one can indeed distinguish a "vortex 
glass" phase. Unfortunately, because of the ever slower 
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timescales over which the spatial structure of the system 
evolves, a precise characterization of this phase using sim- 
ulations of the CGLE is currently impossible. 

To obtain further insight into the problem, we reduced 
the CGLE to a set of ordinary differential equations de- 
scribing the motion of the cores rj (t) and the phases 4>j (t) 
of individual spirals. Only pair interactions are taken into 
account, making use of the results of 0@. The origi- 
nal partial differential equation is thus replaced by a set 
of 3N first-order equations governing the core positions 
Tj(t) and the spiral phases <t>j(t)- The phases determine 
the configuration of the shock lines, which, in turn, yield 
the velocity of the cores. One clear advantage of this ap- 
proach, in addition to virtually suppress physical space, 
is that it allows to bypass the ultraslow timescale related 
to the vanishing value of the wavenumber k as c — * 0. 

The problem of the interaction between two oppositely- 
charged spirals is equivalent to that of one spiral and 
a straight shock line, i.e. a half-plane with the no-flux 
boundary condition dA/dn — 0. For c -C 1, the equa- 
tions governing the position and phase of the spiral are 



dz o / 

— = 2c nkiim 
dt y 



cB 



, exp(-2|cfe|X) 



^/ir\ck\X 



2ck2% eM-2\ck\X) 



(3) 



^w\ck\X 

where z = x + iy is the (complex) spiral position, 
dz/ dt = v x +ivy, <fi is the phase correction, m is the topo- 
logical charge, B' ~ 0.48, and X is the distance from the 
core (located at (— X, 0)) to the boundary (the y-axis is 
chosen to be on the shock line) [ ^2[ . Eqs. ([|) imply the 
asymptotic repulsion of the spiral from the boundary, in 
agreement with simulations JlT[ | . This prevents spiral an- 
nihilation in the vortex liquid. 

The equations in the oscillatory range c > c* read: 



C x v x + mC v v v = 



-kyl — k 2 exp(— pX) 



5^Hp~X 



X' 



(4) 



Im(doQ) 



Im 



-ky/l - k 2 exp(-pX) 

8y/2^p~X 



x-^ 



where the complex constants C x , y , Co, C'iq are obtained 
numerically and the real parameters p, fi, S are derived 
from the linearized CGLE (T^Q. Eqs. (@) describe 
stable bound states of spirals moving along the plane 
boundary. However, the problem is complicated by the 
fact that the position of the shock line (where the waves 
emitted from the cores collide) depends on the relative 
phase of the spirals. From the condition that the total 
spiral phases <&j = ip( r ) +4>j ~ — l^k + ^j are equal, one 
finds that the distance to the shock is given by: 

X = \v k -v ] \/2-( ( p k -^)/2\k\ (5) 

Substituting (|5|) into Eqs. (^) one finds that the symmet- 
ric bound states are unstable in the oscillatory range [O . 
Similar phenomena occurs for likely-charged spirals. 



Eqs. (§), (|) and condition © were used to investi- 
gate the dynamics of many-spiral states. In this case, one 
considers each spiral in the local coordinate system asso- 
ciated with every other spiral and sum up their contribu- 
tions (details will be published elsewhere). Substituting 
r — > 2c\k\r,t — ► A\/2nc 3 k 2 and introducing a rescaled 
"phase length" £ = ccf>, one obtains from Eqs.(0): 



^ = £(-cB' + »n fe ) Zk 



exp(~X k 



\zk - Zj\ 



'Xh 



exp(-X fc ) 
2cJX~l, 



(G) 



where X k — \z k ~ Zj\/2 + Q — ( k - A similar rescaling was 
performed in the oscillatory case. 

We investigated numerically and analytically the re- 
duced equations which, in fact, are rather different in the 
oscillatory range pT| , p^ | and in the monotonic range [jl2| . 
Overall quantitative agreement with the full CGLE dy- 
namics has been found, which constitutes maybe the first 
success for the ideas underlying the concept of "defect- 
mediated turbulence" put forward in the 80's [p|,^6[. 

In particular, we found, in the monotonic range, clear 
"liquid" behavior accompanied by normal diffusion for 
all c values. Each spiral core moves chaotically and the 
velocity distribution is also very close to a Maxwellian. 
The reason underlying the chaotic behavior of this vortex 
liquid is that Eqs. @ do not obey a variational principle 
for any c value, due to the nontrivial form of the pair- 
wise interaction. However, an ensemble of likely-charged 
spirals tend to form a stable hexagonal "Wigner" crystal 
due to mutual repulsion. With periodic boundary con- 
ditions (used here), one has an equal mixture of positive 
and negative spirals, and it can be proved that a square 
lattice of spirals with alternating charges is unstable in 
a large system with respect to long-wavelength pertur- 
bations. This explains the short-range crystalline order 
and local spatial charge separation observed (Fig. |l|c). 
The faithfulness of the reduced equations to the CGLE 
is also testified by the quantitative agreement: found for 
the value of the diffusion constant D. for c = 0.6, the 
domain size L = 512 and vortex number N = 46 the dif- 
fusion D obtained from CGLE is D w 0.0033 and from 
Eqs. is D w 0.0036. The reduced equations do allow 
for an extensive study of the asymptotic (i.e. large-size) 
properties of the vortex liquid. In particular, we find that 
D ex exp(— 2c\k\j yfp) at fixed c value (Fig. |f]f). Taking 
into account that the typical inter-spiral distance varies 
like 1/y/p, this scaling law is in agreement with the ex- 
pression for the spiral velocity as given by Eq. ([|) . 

In the oscillatory range also, the reduced equations 
faithfully reproduce the phenomenology of symmetry 
breaking, intermittent activity at low vortex densities, 
and "vortex glass" formation observed in the CGLE 
(Fig. ||c;,f). The reduced equations provide an interesting 
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framework to discuss the mechanisms leading to inter- 
mittent dynamics in the vortex glass containing liquid 
droplets. For a spiral embedded in a large domain of 
arbitrary shape, Eqs. ^ generally do not have a sta- 
tionary solution Vj = <j)j = for all j (because the com- 
plex constants C x>y and Cio,o are typically not propor- 
tional). As a result, it is impossible to satisfy the condi- 
tions Vj = and <pj — simultaneously. The drift of the 
spirals changes the phases and, therefore, the shape of 
their domains. Eventually, when the phase difference be- 
comes large enough, the shock line boundaries rearrange 
and trigger rapid dynamics within the liquid droplets. 
Since the velocity of the cores is an exponential function 
of the inter-spiral separation, the typical timescale be- 
tween such rapid rearrangement events is exponentially 
large in the spiral size. Thus, in very large systems, where 
the domains occupied by spirals may be comparable with 
system size, the dynamics will be very slow indeed. 

We have found that the transition from vortex liquid 
to vortex glass possesses some features of a first order 
phase transition. Additional numerical studies will be 
needed to precise the structure of corresponding phase 
diagram. Returning to the question of the possibility of 
"true" glassy dynamics |l7j] (as testified, e.g., by some 
aging phenomena), several comments are now in order. 
Preliminary numerical results of the reduced equations in 
very large systems at low vortex density indicate that the 
intermittent space-time activity usually reaches a com- 
plex but self-averaging asymptotic regime, in contrast 
with spin-glass-type behavior Q and in closer agree- 
ment with a structural glass, which can be viewed some- 
times as a very viscous fluid. On the other hand, these 
regimes, as far as they can be studied for the full CGLE, 
do exhibit a slow, possibly aging-like, decay of the num- 
ber of vortices in the system (Fig. ^a). This leaves 
the door open for an actual vortex glass in this fully- 
deterministic, noiseless, disorder-free context — a chal- 
lenge to models of glassy behavior in statistical physics. 

Our work has shown that the multi-spiral "frozen" 
states of the CGLE actually evolve on very long 
timescales. These ultra-slow regimes are relevant, we 
believe, to the many related models and experimental 
situations. A vortex liquid and a vortex glass can be dis- 
tinguished depending on the effective interaction law be- 
tween spirals. These "phases" are dynamical, and are re- 
markably well accounted for, in the case of the CGLE, by 
a reduction of the full problem to a finite set of ordinary 
differential equations governing the vortex cores. Given 
that the typical timescales of vortex motion are several 
orders of magnitude larger than the period of the basic 
oscillations, any experimental observation will require a 
very good control of the apparatus over long times. This 
work was supported by the U.S. Department of Energy 
under contract W-31-109-ENG-38. 
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